Stable topological superconductivity in a family of two-dimensional fermion models 
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We show that a large class of two-dimensional spinless fermion models exhibit topological su- 
perconducting phases characterized by a non-zero Chern number. More specifically, we consider a 
generic one-band Hamiltonian of spinless fermions that is invariant under both time-reversal, T, and 
a group of rotations and reflections, G, which is either the dihedral point-symmetry group of an 
underlying lattice, G = D„, or the orthogonal group of rotations in continuum, G = 0(2). Pairing 
symmetries are classified according to the irreducible representations of T 8 G. We prove a theo- 
rem that for any two-dimensional representation of this group, a time-reversal symmetry breaking 
paired state is energetically favorable. This implies that the ground state of any spinless fermion 
Hamiltonian in continuum or on a square lattice with a singly-connected Fermi surface is always 
a topological superconductor in the presence of attraction in at least one channel. Motivated by 
this discovery, we examine phase diagrams of two specific lattice models with nearest-neighbor hop- 
ping and attraction on a square lattice and a triangular lattice. In accordance with the general 
theorem, the former model exhibits only a topological {p + jp)-wave state, while the latter shows a 
doping-tuned quantum phase transition from such state to a non-topological, but still exotic /-wave 
superconductor. 

PACS numbers: 74.20.Rp, 05.30.Pr, 03.67.Pp, 71.10.Fd 



I. INTRODUCTION 



The field of quantum condensed matter has provided 
us with quantum many-body states that are nothing 
short of amazing. Among the most remarkable are phases 
associated with a new paradigm^ of topological order. 
These topological phases have a number of fascinating 
and technologically useful properties, such as quantized 
Hall conductance and non-Abelian quasi-particles. How- 
ever, the precise conditions for a quantum topological 
phase to arise from a microscopic Hamiltonian are un- 
known. The few known examples of topological order 
currently include the quantum Hall states^, lattice spin 
models due to Kitaev^, lattice versions of the quantized 
Hall effect^, and related models of topological insula- 
tors^ . While the physical systems that may host the non- 
trivial topological phases are quite different, their the- 
oretical descriptions often formally reduce to that of a 
topological superconductor (SC). E. g., the Moore-Read 
even-denominator fractional quantum Hall wave- function 
is equivalent to the mean-field BCS state of a spinless 
(p-l-zp)-wave SC^. Topological insulators and SCs too can 
be related and classified on an equal footing^-^, by notic- 
ing that any SC is an insulator of its Bogoliubov excita- 
tions, whose "band structure" is governed by Bogoliubov- 
de Gennes Hamiltonian at the mean-field level. There- 
fore, understanding topological superconductivity is an 
important issue both due to its many connections to a va- 
riety of seemingly unrelated topological phases and also 
in its own right, e. g., in relation to the recent experimen- 
tal observation of an exotic paired state in Sr2RuO^>i^ 
and proposals for realization of p-wave superfluids in cold 
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Topological SCs, most notably (p + ip) models, have 



been considered in the theoretical literature in great de- 
tail. However, the starting point of all theoretical mod- 
els has been a quadratic mean-field Hamiltonian, with 
a predetermined topological order parameter of interest, 
or equivalently a reduced BCS Hamiltonian with exotic 
interactions that are difficult to imagine being realized 
in the laboratory. Such models are capable of answering 
some key questions related to the properties of a given 
topological phase, but they do not provide much guid- 
ance in the search of Hamiltonians that would host those 
phases. In other words, these models are sufficient to 
produce nontrivial topological order by design, but do 
not shed light on the minimal necessary conditions for 
the emergence of topological order. 

In this paper we prove a general theorem that allows 
us to construct a large family of lattice models that give 
rise to topological superconducting states. We show that 
contrary to a common perception, the nontrivial topo- 
logical phases do not necessarily arise from exotic Hamil- 
tonians, but instead appear naturally within a range of 
simple models of spinless (or spin-polarized) fermions 
with physically reasonable interactions. Our theorem 
is based on examining the BCS free energy of possible 
paired states which is known to be asymptotically exact 
for weak coupling since BCS instability is an infinites- 
imal instability and the use of the Jensen's inequality, 
which ensures that topological phases are often selected 
naturally by energetics. The paper is organized as fol- 
lows: In Section |TT1 we introduce a general microscopic 
Hamiltonian describing spinless fermions in two dimen- 
sional space and present the BCS mean field treatment 
of superconductivity in this model. The topological clas- 
sification of 2D superconductors is reviewed in Section 
IIIBI and the main energetics argument indicating that 
the topological paired states are energetically favorable 
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is proven in Section III CI In Section IIIII we quantita- 
tively study phase diagrams of two specific lattice fermion 
Hamiltonians on a square lattice and triangular lattice 
with nearest-neighbor hoppings and interactions. The 
ground state of the square-lattice model is proven to be 
a topological {p + ip)-wave SC at arbitrary filling. The 
triangular lattice model gives rise to a (p -I- ip)-wave su- 
perconducting state guaranteed at low filling, but shows 
a first-order phase transition into a non-topological /- 
wave SC at intermediate fillings and another transition 
within the /-wave superconducting dome from a gapless 
to a fully gaped superconductor. 



II. SPINLESS FERMION 
SUPERCONDUCTIVITY IN TWO DIMENSIONS 

We start our general discussion with the following 
single-band Hamiltonian for spinless fermions 



r 4 - - 

Jkk' ,qCi(._|_qC_j^C_k' Ck'+q i 



k6BZ 



q/2,k,k'eBZ 



(1) 



where cj^./ Ck are the fermion creation/annihilation op- 
erators corresponding to momentum k, "BZ" stands for 
"Brillouin zone," = Ck^A* with ek being the dispersion 
relation of the fermions and /i the chemical potential, and 
/k,k',q describes an interaction, which is assumed to have 
an attractive channel. 

We assume that Hamiltonian ^ arises from a real- 
space lattice or continuum model and is invariant with 
respect to the underlying spatial symmetry group, which 
we denote as G, and the time-reversal group, T. We 
note that in two dimensions (2D) the range of possible 
spatial groups, G, is limited to the following dihedral 
point-symmetry groups: Z^i, D2, -D3, D^, and Dq in the 
case of a lattice or orthogonal group of rotations 0(2) = 
Dao in continuum. We recall that the group Dn includes 
^ — rotations and in-plane reflections with respect to n 
axes. The superconducting order parameter is classified 
according to the irreducible representations of the full 
group T(8G. Since, T = Z2, Z2(^L>i = £'2 andZ2®£i3 = 
Dq, we can confine ourselves to studying representations 
of D2, D4, Dq, and 0(2), which exhaust all physically 
relevant possibilities. 



A. BCS Mean Field Theory 

Now, we define the superconducting order parameter 

as 



Ak = J ./k,k'(c-k'Ck'), 

k'eBZ 



(2) 



where /k,k' = (/k,k',q=o - /k,-k',q=o)/2 is the antisym- 
metrized BCS coupling strength. Using a Hubbard- 



Stratonovich decoupling in Eq. ([T]) with q = and ig- 
noring superconducting fluctuations, we arrive at 



^^"= / (ek4^k+iAkctcU + iA 



keBZ 



kC-kCk 



(3) 



k, k'eBZ 



with /i7k' being the matrix inverse of /k.k'- By integrat- 
ing out the fermions we find the BCS free energy func- 
tional expressed in terms of A. It contains two parts, 
T[Ak]=Ti+TH, with 



.F/[Ak] =-T f 
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Tu [Ak] = -2 
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Ak/k:k'Ak'. 



e^ + |Ak| 



(4) 
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B. Topological Classification of Two-dimensional 
Superconductors 

To describe the topological properties of a SC state, 
we introduce the topological index (the Chern number) 
as follows^ 



C 



An 



m • dkm X dkm. 



(6) 
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where m = (mi, 7712,7713) — (ReAk, — ImAk, ^k)/£'k and 
i?k = Vsk + |Akp. This topological index classifies all 
maps from to representing the unit vector m(k) 
into equivalent homotopy classes. We will call a SC state 
topological, if C ^ 0. 

The Chern number is equal to the sum of the winding 
numbers, C = X^o-^o-' which can be defined for each 
segment of the Fermi surface (FS), Va, as follows: 



27rW^ = * VkV3k • dk 



(7) 



where ip-^ is the complex phase of Ak. Note that even 
though we assume a single-band picture, a general sit- 
uation is allowed where the FS is formed by one or 
more disconnected components, FS — Va with a = 

1,2.. ..,77. 

To prove the relation between C and Wo-'s, we separate 
the closed Brfllouin zone, d{BZ) = d {S^ x S^) = 0, into 
an "electron" region, Ebz — {k. G BZ : 7773 (k) > 0} and a 
"hole" region, Hbz = {k e BZ : 7773 (k) < 0}. The Fermi 
surface is a directed boundary of these regions, FS = 
J2cy = dEBz = —dHsz- One can show that 
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Eq. jl]) and the Stoke's theoremlS yield C = J2aW<y- 

If Wa = for all tr, the complex phase of the pairing 
order parameter can be gauged away via a non-singular 
redefinition of the fermion fields and corresponds to a 
topologically trivial state. This however is impossible if 
at least one winding number is non-zero. We will call such 
states time-reversal-symmetry breaking (TRSB) states. 
The class of TRSB superconductors is larger than and 
includes that of closely related topological SCs. If there 
is just one singly-connected FS, the two types of states 
are equivalent. 

C. General Theorem of the Stability of TRSB SC 
States 

Now we examine the stability of TRSB SCs. The order 
parameter in a certain channel corresponding to a dr- 
dimensional irreducible representation, T, of the group 
T (8) G can be written as a linear combination of real 
eigenfunctions of F, ^^(k) (with a = 1, . . . , c?r) 

dr 

Ak = 5]Aa0^(k). (9) 

a=l 

In two dimensions, the number of irreducible represen- 
tations to be considered is highly constrained and in- 
cludes only ID and 2D real representations. In particu- 
lar: (i) For a system with a four-fold rotational symme- 
try (e. (?., arising from a square lattice), the correspond- 
ing point group, D4, has only one space-inversion-odd 
irreducible representation, E, which is two-dimensional; 
(ii) With a six- fold rotational symmetry (e. g., due to 
a triangular or hexagonal lattice), there exist three irre- 
ducible representations of Dq odd under space inversion: 
A 2D representation, Ei (corresponding to a p-wave pair- 
ing) and two ID representations, Bi and B2 (correspond- 
ing to two types of /-wave pairing), (iii) The continuum 
group, 0(2), has an infinite set of 2D real representations, 
classified by odd orbital momenta, / = 1,3,5,.... 

We now consider a pairing channel corresponding to a 
2D representation of T ® G. There are two real eigen- 
functions for this representation: (pi (k) and (j)2 (k) . If the 
order parameter is proportional to either of them, it is 
real and corresponds to a topologically trivial state with 
zero winding number. We prove below that such a state 
is always unstable. The invariance of the Hamiltonian 
under T (g) G ensures ^Non-top = ^ ['/'i(k)] = [4'2(k.)] 
(e. g., Px- and Py-states have the same energies in con- 
tinuum) . 

Let us show that one can always construct a new TRSB 
state with 

0TRSB(k) = [01 (k) -I- i02(k)] 

that has a lower free energy than .^Non-top- One can 
see from Eq. (O that [4>t[isb0^)] = ^11 [4>iO^)] = 
Tu [02(k)] because 0|RSB(k) = 0?(k)/2 + 0i(k)/2,. To 



handle the less trivial "quasiparticle part" of the free 
energy dU we take advantage of the Jensen's inequal- 
ity which states that for any function with fix) < 0, 
f{x/2 + y/2) < f{x)/2 + f{y)/2 for any x ^ y. The 
integrand in Eq. ([4]) for is a concave function of 
X — |Akp and therefore satisfies the Jensen's inequal- 
ity (which after integration over momentum becomes a 
strong inequality for all physically relevant cases). Since 
•^trsbC*) = + '/'2(l^)/2, we have proven that 

TU ^ •^[0iW]+-^[<^2(k)] _^ 

■r [(pTRSBikjJ < ^ = >Non-top- 

(10) 

This inequality (to which we refer to as "theorem" ) repre- 
sents the main result of our work and proves that a TRSB 
phase is always energetically favorable within a 2D repre- 
sentation. This is a strong statement that is completely 
independent of microscopic details, such as hoppings and 
interactions, and relies only symmetry. It leads, in partic- 
ular, to the conclusion that any single-band spinless SC 
(and certain models of spin-polarized SCs) originating 
from a square lattice with singly-connected FS must be a 
(p-l-ip)-paired state. Similarly, any SC arising from spin- 
less fermions in continuum must be of a {2l + l)+i{2l + l)- 
type, which is topologically nontrivial. This includes 
all continuum models with attractive forces and con- 
ceivably some continuum models with weak repulsion 
that may give rise to pairing via Kohn-Luttinger mech- 
anism— li^ii^. Since a large number of lattice fermion 
Hamiltonians at low particle densities reduce to an effec- 
tive single-band continuum model, it means that at least 
in this low-density regime any paired state is guaranteed 
to be topological. 



III. LATTICE MODELS 

To illustrate how our theorem manifests itself in prac- 
tice, we examine specific models within a large class of 
generic tight-binding Hamiltonians on a lattice 

r.r' r (r,r') 

where cj/cr creates/annihilates a fermion on a lattice 
site r. We note that this real-space Hamiltonian re- 
duces to a more general model ([1]) via a lattice Fourier- 
transform. For the sake of concreteness, we focus below 
on the following two models with nearest-neighbor hop- 
pings, ir,r' — tS\r-r'\.i E^^d nearest-neighbor attraction, 
Vry = —gS\r-T'\,i on (i) a simple square lattice and (ii) a 
simple triangular lattice. 



A. Square Lattice 

The square lattice case corresponds to the D4 sym- 
metry group, which has only a 2D representation. The 
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FIG. 1: (Color online) The phase diagram for fermions on a 
square lattice with nearest-neighbor hoppings and attraction 
{g/t=l). The phase boundary separates a normal metal and 
a topological {px + ipy)-waYe SC. The insets display FSs for 
^J, <0 (left) and fj, > (right). 



attractive interaction guarantees that the ground state 
is a SC^^ and the general theorem pU]) guarantees that 
it is topologically non-trivial. To see how this happens 
in the specific model, we define two independent or- 
der parameters on horizontal and vertical links: A„ = 
g{crCr+e„), whcrc n = X OT y and e„ is the correspond- 
ing lattice vector (we use units where the lattice con- 
stant, a = 1). These real-space order parameters are 
related to the momentum-space definition Q via Ak = 
2i Aa4>aO^), with the BCS interaction being /k,k' = 

a—x^y 

—g S 0a(k)(/'Q(k'). Here we defined two eigenfunc- 

ct—x,y 

tions of the above-mentioned 2D representation of Dj^: 
(j)x.yi}i) — sin (k • e-x.y). It is straightforward to calculate 
the BCS free energy given by Eqs. @ and ([5]) for all pos- 
sible order parameters encompassed by the linear com- 
binations Ak = g [Aa:0a;(k) -I- Ay0y(k)], wlth arbitrary 
^x,y S C. We find that a [p + ip)-superconducting state 
with \x — ±iXy is selected at all /i. Fig. [T] summarizes 
the phase diagram of the model on the fJ- — T plane. The 
maximum Tc within the mean-field treatment occurs at 
half-filling. The tails of the particle-hole symmetric phase 
boundary correspond to small "electron" and "hole" den- 
sities, and therefore to continuum limit with the isotropic 
quadratic dispersion, ^k — (fc^ — fcp) /(2m*), the effec- 
tive mass, m* = l/(2to^), and the Fermi momentum, 
kpa = ^|^±4t|/(2i). 

It is useful to consider the continuum limit 



l/i ± At\ /t ^ in more detail, as it gives a valuable in- 
sight into stability of the topological phases. For this 



Eqs. (HI) and ([5]) to derive the Ginzburg-Landau free en- 
ergy (per unit area): 



A 



where ly — m* /{2tt) is the density of states at the FS, 
Tc is the BCS transition temperature, C is the Riemann 

zeta- function, Aq = g^JlX^f^ + |Aj,|^ is the modulus of 



the order parameter, A is the area of the sample, and 
we introduced a symmetry factor, S, as follows [below, 
6'k = tan"^ (ky/kx)] 



d^k IXxM'^) + >^yMW 



2tt 



ksFS 



iA.r + iA,r 



(12) 



The minimal free energy below Tc is given by 
^GL,mi„/-A = -5„,ax [47ri.T2 lu^ (T/Tc)] /[7C(3)]. There- 
fore, the absolute minimum is achieved by maximizing the 
symmetry factor, S. In the continuum limit |k|a — > 0, 
we can approximate the normalized eigenfunctions of 
D4, by (pxO^) = ■\/2cos(0k) and 0y(k) = \/2sin(6'k). 
Hence, the topologically trivial px- and p^-states lead to 



(4< 



1 

FS 



2/3, while the topological states 



Px ± ipy yield 5p^±ip„ = (^|e*'''''| = 1 > 2/3 and 

therefore are selected by energetics. This fact is a special 
case of our general theorem summarized by Eq. (jlOp . 

We note that the mean-field BCS-type model can for- 
mally be considered for the extreme values of the non- 
interacting chemical potential > 4i, which is not as- 
sociated with a non-interacting FS. Hence, mean-field 
paired states in this limit are not topological and corre- 
spond to the strong-pairing (Abelian) {p + ip) -phase con- 
sidered by Read and Green^. While such a mean-field 
BCS model is sensible in the context of the quantized 
Hall state, it may be unphysical for fermion lattice mod- 
els. Indeed, the chemical potential, /i, is renormalized by 
non-BCS interactions or equivalently by superconduct- 
ing fluctuations originating from the terms with q 7^ in 
Eq. ([1]). These strong renormalizations are bound to shift 
fi towards the physical values with a reasonable Fermi 
surface, which in a metal is guaranteed by Luttinger the- 
orem. Hence, it is not clear whether the Abelian (p -\- ip) 
superconducting states may survive beyond mean-field. 
Due to these arguments, we disregard here such case of 
non-topological [p -\- ip)-paired states. 

We now derive Bogoliubov-de Gennes equations from 
the lattice model. These equations are often the starting 
point of discussions on bound states in a vortex cor o^^i^° 
and edge states^. To do so, we first present the fermionic 
mean-field BCS Hamiltonian on a lattice as follows: 



MP 



r> hr 



Ar 



which is a real space version of Eq. Q where Arr' = 
g(crCr') is the order parameter on the bond (rr') and 



purpose, we use standard perturbative expansion's in "rr- 



-M|r-r'|,i ~ fJ-Srr' IS the matrix element of the 



single-particle Hamiltonian. We then follow the stan- 
dard route and introduce Bogoliubov's transform — 
7Ur + commutation relation [7i'^^,7] — 

—Ej. This yields gives the desired BdG equations 



EVr = ^ (-AJr'Ur' - ft-rr'^r') 



(13) 
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In principle, the order parameter Arr' should be deter- 
mined via solving BdG equation self-consistently. How- 
ever, we know that in a homogeneous ground state the 
order parameter has a (p -f ip)-wave pairing symmetry, 
i.e., Ay — ±iAa;. If there are inhomogeneities in the 
system (e.g., vortices, domain walls) the pairing symme- 
try (associated with the relative phase between Aj, and 
Aa: components) is not necessarily p -I- zp. But since this 
pairing symmetry is selected by energetics, we expect 
such deviation to be irrelevant for low energy physics. 
Therefore we can assume that the relation A^ = iiAj, 
holds for general configurations of order parameter at the 
mean-field level. This is equivalent to separation of the 
Cooper pair wave function into parts corresponding to 
the center-of-mass motion and relative motion. 
Now we take the continuum limit of 

where m* is the effective mass and /i = -I- 4f is the 
chemical potential measured from the bottom of the 
band . To treat the off-diagonal part, we formally 
represent the second term in Eq. ([TS]!) as follows 
Er' ^rr'Wr' — Ai;(r), with the gap operator being 

A = ^Arr'e('"'-^'-^-. (14) 

r' 

The order parameter Arr', which "lives" on bonds, 
should be castcd into only site-dependent form as fol- 
lows: 

Arr' = A ^^P(*^r'-r) (15) 

where 0r'-r is the polar angle of r' — r. Then, we expand 
(fn)l to first order in |r' — r| = a and obtain the familiar 
BdG equations in continuum: 

(16) 

Ev[v) = AtM(r) - ^[-iV)v[v) 

where the gap operator A = a{A(r),9a; -I- idy\. An in- 
teresting question to be addressed elsewhere is whether 
fluctuations and in particular deviations of pairing sym- 
metry from p+ip play a role in the topological properties. 

B. Triangular Lattice 

We now address the very interesting case of a simple 
triangular lattice. Here the Dq symmetry group has both 
a 2D representation (p-wave) and two ID representations 
(/-wave). Therefore, non-topological /-wave states are 
allowed. Low "electron" densities correspond to a single 
circular-shaped Fermi surface and must lead to the p+ip- 
wave pairing per the same argument as above. However, 
the spectrum of the model is not particle-hole symmetric 
and at large fillings (with /i > /x* = 2t), the electron 
Fermi surface splits into two hole-like Fermi pockets and 




-6 -3 3 



fl/t 

FIG. 2: (Color online) The phase diagram for spinless 
fermions on a triangular lattice with nearest-neighbor hop- 
pings and attraction {g/t=l). The bottom of the band is 
located at = —6t and the top is at = 3t; fi* = 2t cor- 
responds to a van Hove singularity. Two SC phases, with 
{Px + iPy)- and /-wave symmetries are present. They are sep- 
arated by a first-order phase transition at fia/t ~ 1.057. The 
insets (left to right) are the FSs for fi < ^* , n < fi* , and 
^ > H* and the dashed lines indicate the nodal directions of 
the /-wave SC. 

maps onto an effective continuum model but with two 
fermion species: 

^2h,off = y (Ck^+^k'^+.k + ^-k/j-L_k^-,k) + interactions, 

k 

(17) 

where /i±,k are fermion operators near the two pockets 
labeled by a pseudospin index ct = ± and the spectrum 
is asymptotically given by 

i]^ = k^/2m + a{kl-ik^kl)- Ey, (18) 

with k measured from the corner points of the hexago- 
nal Brillouin zone. Note that under a tt or ±7r/3 rotation, 
the spectrum transforms as ^k ^ C-k and this symmetry 
is preserved if cr — *■ —a. This leads to a pairing analo- 
gous to the s-wave pairing of spin-1/2 fermions, with the 
order parameter of the inter-pocket pairing defined as 
A/i = g J'j^(/i_|__k^-.-k)- However, this is an /-wave pair- 
ing state, because under a 7r/3-rotation, a — > —a and 
A;i(k) changes sign. Since the low-density limit leads to 
a topological phase and the high-density limit leads to an 
/-wave topologically trivial state, there must be a quan- 
tum phase transition in between. The entire phase dia- 
gram can be derived using Eqs. ([4]) and ([5]) and the real- 
space construction as follows: On a triangular lattice we 
can define three order parameters on the nearest neighbor 
bonds corresponding to the three lattice vectors, e„ with 
azimuth angles 2n7r/3 and n = 0, 1, 2: A„ = g{crCj.^e„) ■ 
Two different types of pairing channels are formed by 
these three order parameters: An /-wave channel with 
A„ = A and a p-wave channel with A„ = Ae^^'^™/'^ . 

The resulting phase diagram is shown in Fig. [2l As 
expected, a topological p + ip-wetve SC state with A„ = 
^g±2irin/3 jg stabilized at low fillings, while an /-wave 
state with A„ = A is favored at high densities. These 
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phases are separated by a first-order transition. As 
shown in Fig. ^ the van Hove singularity fj, = fj,* gives 
rise to a maximal Tc and is located inside the /-wave 
superconducting dome. This point represents another 
type of a quantum transition that separates two qualita- 
tively different topologically trivial paired states: (1) For 
fi < fj,* , there is just one electron- type Fermi pocket 
that is cut by the nodes of the /-wave gap in the di- 
rections, o'^l^ — mvr/S + 7r/6. This gives rise to gap- 
less quasiparticles. (2) For fi > fi* , no FS can be cut 
and the nodal quasiparticles disappear. The phase be- 
comes fully gapped and eventually crosses over to the 
two-specie continuum model (|17|) . Experimentally, the 
two types of /-wave phases can be distinguished by dif- 
ferent T-dependence of the heat capacity. 

We also present the Bogoliubov-de Gennes equations 
for the /-wave pairing state in high-density limit fi 3t. 
Their derivation goes along the same lines as that given 
in Section FlII Al for v^+ivy pairing SC. However, in the /- 
wave case, the momentum space order parameter is given 
by Ak = A(sink • ei -I- sink • 62 + sink • 63). Therefore, 
the order parameter reads: 

Arr' = A (^^-^^ C0s(36lr'-r). 

Since 61-1-62 + 63 = 0, the leading term in the expansion 
is ~ a'^. With some algebra one can show that the gap 
operator is 

^^^j2i9nAdn,{dn,^{r)}}} (19) 

with 9„ = V • e„ and the BdG equation takes the form 
of Eq. (Uni). 



IV. DISCUSSION AND CONCLUSIONS 



In conclusion, we discover that topological supercon- 
ducting phases breaking time -reversal symmetry emerge 
naturally within a large class of spinlcss fcrmion models. 
The technique we apply here has a close relation to BCS 
mean field theory of a spin-triplet superfluid ^He^'^''^'^, 
which concluded that the B-phase with isotropic gap is 
stabilized compared to anisotropic A-phase^^. However, 
we have shown that similar conclusion can be generalized 
to any band structures, filling factors, and interactions, 
as long as the system satisfies proper (discrete) rotational 
group symmetries. More importantly, our proof is insen- 
sitive to the existence of the "nodes" . In continuum, it 
has been argued that a Px state is unstable against the 
Px + ipy pairing state, because the former has nodes thus 
having smaller condensation energy. However, the sta- 
bility of a nodeless Px state, which could exists in lattice 
models, was unclear before this our work. 

We should also emphasize that although the discus- 
sions above focus on spinless fermions, all the conclusions 
can be generalized to the triplet pairing channels of spin- 
1/2 fermions, because these pairing channels also corre- 
spond to the space-inversion odd representations of the 
symmetry group. In addition, we note that any pairing 
state that spontaneously breaks a lattice rotational sym- 
metry must have at least one degenerate state for both 
spinless and spin-1/2 fermions. Our theorem indicates 
that these type of states must have a complex pairing 
order parameter to be energetically stable. 
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